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Table 2 — Continued 
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Table 2 — Continued 
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1. Introduction 

Mira-typc variable stars are large, cool stars whose visual light variations exceed 2.5 
magnitudes over periods from 150 days to ~500 days. The light curves of Mira variables 
depend on the surface temperature, radius, and opacity, all which vary as the star pulsates. 
These pulsations extend the atmosphere beyond that of the hydrostatic equilibrium 
configuration and enhance mass loss in these stars (Maciel 1977; Willson & Hill 1979; 
Bertschinger & Chevalier 1985; Bowen 1988; Fleischer et al. 1992, 1995). As a result, Mira 
variables are an important component in seeding the interstellar medium with C, N, and O. 

These stars are located on the asymptotic giant branch, a transitional phase in stellar 
evolution. Photometric and spectroscopic measurements of their light curves provide a 
means to probe the stellar atmospheres and underlying pulsation mechanisms occurring 
during this stellar phase. In the near infrared, the spectra of Mira stars are dominated by 
the TiO 7 system, the VO 7 system and ZrO (Wing 1967). The TiO features are thought to 
be produced in a layer somewhat far from the photosphere (Gillet 1988). HanifF et al. (1992) 
present optical aperture synthetic images of the photosphere of Ceti at 6500 A, 7007 A, 
and within a TiO bandhead at 7099 A, with the star phase ~0.94. They find asymmetry in 
the images, with the TiO image one and a half times larger than the photospheric images. 
Also, narrowband speckle interferometric measurements taken in the TiO 7120 A bandhead 
and outside at 7400 A by Labeyrie et al. (1977) shows that the diameters of R Leo and o Get 
are twice as large in the TiO feature than outside of it. This demonstrates that a model 
atmosphere, based on the spectra observed over a TiO bandpass, provides parameters such 
as Tgff and log(5f) in an atmospheric layer far from the photosphere. 

Joy (1926) presents a comprehensive phase dependent spectroscopic (35 A/mm) study 
of a prototype Mira variable, Ceti. Analyzing 131 spectra taken over a ten year period, 
Joy describes several important characteristics. Briefly, the spectra of Ceti show that 
TiO bands vary with magnitude, hydrogen emission lines appear with greatest intensity 
at or shortly after maximum visual brightness (phase zero), ionized iron emission lines 
are observed at maximum, and the temperature is estimated to vary from 1800 K to 
2300 K. Absorption lines (including iron, vanadium, chromium, manganese, calcium, and 
magnesium) were used to measure a variation in radial velocity. The maximum positive 
velocity occurs at phase zero, and greatest blueshift at minimum light. Later, Joy (1954) 
took 88 spectrograms of Ceti at a higher spectral resolution (typically 10.3 A/mm) and 
confirmed that maximum velocity of recession occurs soon after visual maximum. He 
attributed these results to a pulsational mechanism, and suggested the possibility of shocks. 

Perhaps one of the most interesting characteristics of Mira spectra is the strong 
hydrogen Balmer line emission that is seen throughout much of the pulsation cycle. As 
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Pickering (1887) first noticed in spectra of o Ceti, and later described in detail by Joy (1926, 
1947, 1954), the hydrogen Balmcr emission line flux in relation to the nearby photospheric 
(i.e., pseudocontinuum) flux is unique in the oxygen-rich (M-type) Mira spectra: Balmer 
Ha emission is typically weaker than H/3 which in turn is weaker than H7 near peak visual 
brightness. is seen as the strongest Balmer emission line at phase zero. Lines higher 
in the series (i.e., towards shorter wavelengths), are weaker. This Balmer increment (i.e., 
FHa < Fh/3 < Fff-y < Fhs) is just opposite of what would be expected. Ha having the 
largest oscillator strength, should be stronger than H/3 and the higher order Balmer lines 
should be weaker down the line (i.e., one should see a Balmer line decrement), assuming 
these lines all form in the same region of the atmosphere (i.e., similar T and P). Meanwhile, 
in S-type and carbon-star (N-type) Miras, the strength of the Balmer lines approximately 
follow their expected respective oscillator strengths (Merrill 1940). 

For years, this Balmer-hne increment in the M-type Miras has been attributed to 
TiO absorption which hide Hct, H/5, and H7 fluxes (Merrill 1940; Joy 1947; Gillet 1988), 
although there is some debate of the extent that this or other molecular absorption has 
on the Ha line (Gillet, Maurice, & Baade 1983; Gillet et al. 1985). Recently another 
explanation has been given for this Balmer increment: NLTE radiative transfer calculations 
of hydrodynamic models representative of Mira variables (Bowen 1988) suggest that the 
Balmer increment results from radiative transfer effects in the hydrogen lines themselves 
when formed in a shocked atmosphere (Luttermoser & Bowen 1992; Luttermoser, Bowen, 
& Willson 2000, in preparation). In these calculations. Ha, having the highest optical 
depth, forms just in front of the innermost shock. H/3 then forms a little deeper, due to 
its lower optical depth, and H7 deeper still. The optical depth of H(5 causes it to arise 
from the hottest part of the shock, hence extends higher above the continuum then the 
longward Balmer lines. Then as one goes to higher-order lines in the scries, the opacity in 
these lines is not high enough in the shock for these lines to form there — we sec through 
the shock at these transitions. As such, this increment may be giving us information of 
the shock thickness for Miras and may indicate that the shock structure of the S-typc and 
N-type Miras are fundamentally different than the M-type Miras. Future NLTE radiative 
hydrodynamic models of Miras with different C/0 ratios are needed to see if this is the 
case. It is likely that a combination of both processes (i.e., TiO absorption and NLTE 
radiative transfer effects) are responsible for this Balmer line increment. 

Another striking feature is the weakness (and often absence) of the He hue (3970.074 A) 
in the Balmer series near maximum visual brightness (see Gillet 1988). Merrill (1940) noted 
this and suggested this weakness a result of the interaction between the He transition and 
the Ga II H line (3968.470 A) wing. Gastelaz & Luttermoser (1997, hereafter CL97) concur 
with this suggestion. Briefly, the He photons may be scattered by the Ga II H line out to 
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IR wavelengths via the 8662 A hne. Of the three hnes in the Ca II IR triplet, the 8498 A 
and 8542 A lines share the same upper level of the Ca II K line, whereas the 8662 A line 
shares the same upper level as the Ca II H line. As such, if He photons arc being scattered 
by Ca II H, the 8662 A will show variations independent of the other two Ca II IR triplet 
lines and can be tested by monitoring the strength of the absorption of this 8662 A line as 
compared to the 8498 A and 8542 A lines as a function of phase in the Mira stars. 

Ca II IR triplet observations of Miras have been carried out by Contadakis & Solf 
(1981) and Gillct ct al. (1985). Unfortunately, Contadakis & Solf (1981) only made one 
observation of the A8662 line in their monitoring program of S-type Miras. They observed 
the A8498 and A8542 in emission in many of the stars in their sample near phase and 
as absorption lines at other phases. So, we cannot use their data to test our proposed 
hypothesis. Gillet et al. (1985) observed P Cygni profiles in the Ca II IR triplet hnes near 
phase in the hot Mira variable, S Car. Once again, only one spectrum of the 8662 A 
line was obtained in their sample of Ca II spectra. CL97 present a set of phase dependent 
spectra taken especially to address this problem of He photon scattering by Ca II. Spectra 
of seven Mira variables were taken at different phases and suggest a possible anticorrelation 
between Ha emission and Ca II A8662 absorption. Assuming that the He line strength 
variations are in phase with Ha, then the apparent anticorrelation between the strength of 
the Ha emission line and the strength of the Ca II A8662 line suggests that a fluorescence 
is taking place in the Ca II A8662 line with He serving as the pump through the Ca II H 
line. In this paper we present phase dependent spectra of twenty Mira variables to further 
explore the anticorrelation of Ca II A8662 absorption with Ha emission. 

2. Observations 

The twenty stars for which spectra are presented, their equatorial coordinates, mean 
spectral types, and visual maxima and minima, taken from the SIMBAD database are 
presented in Table 1. Also in Table 1 are new ephemerides of these Mira stars. We 
calculated the ephemrides from AAVSO light curves measured within the previous eight 
years. The curves were fit with a linear combination of sine and cosine functions from which 
a new period and JDO were determined for each star. 

Spectra of Mira variable stars taken between February 1997 and November 1998 using 
a low resolution spectrograph. Since we are interested in only monitoring absorption and 
emission line strengths for this program, low-resolution spectra are all that is required. The 

spectrograph was used at both the Southeastern Association for Research in Astronomy 
(SARA) 0.9-m telescope at Kitt Peak, and Appalachian State University's Dark Sky 
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Observatory (DSO) 0.45-m telescope located near Boone, NC. A converter lens was used at 
both sites to convert the respective telescope f-ratio to about f/11 for the spectrograph. 

The spectrograph was configured with a 600 g/mm grating. The sht width was 100 
fj,m. At the SARA 0.9-m telescope, the slit width is 3 arcseconds, and at the DSO 0.45-m 
telescope, the slit width is 4 arcseconds. At both sites, the slit was parallel to the hour 
angle. A cooled 768 x 512 CCD camera with 9 //m square pixels was used to record the 
spectra. The linear dispersion is 1.08 A per pixel, covering 768 pixels, or 829 A on the CCD. 
The spectral resolution was measured to be 2.4 A. By rotating the grating up to four times 
per star, spectra were taken from about 6200 A to 9000 A and include Ha, TiO, VO, and 
the Ca II infrared triplet lines. Integration times were adjusted to achieve a signal-to-noise 
greater than 100 for most spectra, except for the May 1997 spectra of W Her, T Her, and 
W Lyr, and the May 1998 spectrum of T Her. 

Table 2 gives the log of observations for each star, which includes dates of observation, 
phase of the variable, approximate visual magnitude, integration time, and observing site. 
The phases listed in Table 2 were determined from the ephemerides given in Table 1, and 
refer to the visual phases, with phase zero corresponding to maximum visual brightness. 
The visual magnitudes have been obtained from curve fits to the light curves from the 
AAVSO database. 

Dark frames and fiat frames were taken for fiat fielding purposes. Spectra of neon 
lamp emission were taken simultaneously with the stellar spectra and used for wavelength 
cahbration. We fiat-fielded the images and extracted the spectra using MIRA software. 

The extracted spectra were wavelength calibrated using the spectrum of neon superimposed 
on the CCD frame with the stellar spectrum (Crowe, Heaton, & Castelaz 1996). 

3. Results 

The spectra for the Mira variable stars are shown in Figure 1. The wavelengths of 
TiO, VO, ZrO, the Ca II IR triplet. Ha, and terrestrial oxygen are marked above each 
set of spectra. Their wavelengths were identified in CL97. The ZrO and VO absorption 
overlap at 6574 A and 6578 A, and VO and TiO absorb at 7865 A and 7861 A, respectively. 
Due to the low dispersion of our spectra, these features are blended. The appearance of a 
relatively narrow feature near 8230 A is seen in some of the spectra - an o Cct spectrum 
taken 1 November 1997 (phase 0.80), and spectra taken 26-28 May 1997 of R Leo (0.49), 
V CVn (0.65), R Boo (0.11), R Ser (0.26), RS Her (0.75), and W Lyr (0.94). This feature 
is due to terrestrial H2O at 8227 A (Turnshek et al. 1985) and is an effective measure of the 
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relative humidity in the air. 

3.1. Radiative Transfer in the Ca II Ion 

The transfer of radiation in Ca II ion is very complicated. Beside He photons affecting 
the level populations of the 3d ^Pi/2 state through the Ca II H line, the hydrogen Lyman-a 
hne hes just shortward of the Ca II ionization edge of the 3d ^03/2 (at 1218.1 A) and 
the 3d ^D5/2 (at 1219.0 A) states. These two states are the lower levels of the Ca II IR 
triplet lines and arc mctastablc. If Lyman-a is a strong emission feature, and the Ca II 
continuous opacity forms in a region of the atmosphere where the Lyman-a line is not yet 
in detailed balance, then Lyman-« photons may influence the Ca II IR triplet lines as well. 
Note that no observations have yet been made in the far-UV for Miras. 

Line center of Lyman-a lies 2.4 A shortward of the ^03/2 ionization edge and 3.3 A 
shortward of the ^05/2 edge. So, much of the Lyman-a emission profile can ionize Ca II 
out of the mctastable state. The lower level of the Ca II 8662 A line transition is the 
^D3/2 state whose ionization edge lies slightly closer to Lyman-a than the ^05/2 edge. 
The following question arises: Will Lyman-a photons affect the level densities of the two 
3d states differently? To answer this question, we ran a few atmospheric models with 
an arbitrarily located 10,000 K shock through the LTE stellar atmosphere code ATLAS 
(Kurucz 1970; Brown et al. 1989). Although it is hkely that NLTE effects will dominate 
the level and ion densities in the atomic and molecular species (Luttermoser & Bowen 
1992; Luttermoser, Bowen, & Willson 2000, in preparation), these LTE runs are performed 
to merely determine the variation of the Ca II continuous opacity across the Lyman-a 
profile. It also should be pointed out that in regions of the atmosphere where LTE no longer 
applies, the assumption of radiative and hydrostatic equilibria are no longer valid either in 
these pulsating giant stars. In fact, a Mira star has numerous shocks propagating through 
its atmosphere at any given time as has been shown by Bowen (1988) and more recently by 
Hofner, et al. (1998) and Loidl, et al. (1999). Willson (2000) gives a very detailed review of 
all the dynamic modeling that has been performed on these pulsating stars and discusses 
the problems of carrying out NLTE radiative transfer in such a dynamic atmosphere. 

We sampled atmospheric depths in front of the shock, in the shock, behind the shock, 
and deep in the photosphere where the continuum reaches optical depth unity in this region 
of the spectrum. We found that the continuous opacity from the Ca II 3d ionization 
remained constant (from both J sub-levels) to within 0.15% from the location of the 
edges through 1210 A, which should include most of the Lyman-a emission profile. The 
fact that the continuous opacity from Ca II remains constant across the Lyman-a profile 
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indicates that Lyman-o; will not preferentially affect the number density in the ^03/2 level 
as compared to the ^05/2 level — photoionizations of Ca II due to Lyman-o; photons will 
affect the strengths of three Ca II IR lines in a similar fashion. Therefore, any variation in 
the 8662 A line that is not seen in the other lines must result from some process other than 
Lyman-o; photoionizations. 

3.2. The Ua Emission Line and Ca II IR Triplet 

Wc arc interested in the strength of the Ca II IR triplet compared to the strength of 
Ho; as a function of phase, since we are assuming that variations in He will mimic variations 
in Ha. The Ca II IR triplet lines are not strong, as expected for stars later than MO (Zhou 
1991). A total of twenty-seven spectra of fifteen stars in our sample span the wavelength 
range from Ho; through the Ca II IR triplet. The remaining spectra are missing either the 
Ho; or the Ca II IR triplet regions of the spectrum because the spectrograph grating was 
not rotated sufficiently during observation to cover those parts of the the spectrum. 

The observations which show obvious Ho emission features include U Ori (phase 0.29), 
R Leo (0.34), V CVn (0.12), R CVn (0.09 and 0.17), R Ser (0.26), W Lyr (0.94), and x Cyg 
(0.04). At the same phases, the Ca II A8662 is seen in emission in R Leo, R Ser, W Lyr, and 
X Cyg, whereas the A8498 and A8542 fines stay in absorption or are not apparent. Ca II 
A8662 is not seen in emission in any other spectra that we took, only in those that show 
Ho; emission. The Ca II IR triplet can be seen in absorption in the remaining twenty-three 
spectra (although in some cases weakly), except for U Ori (phase 0.01), S CrB (0.48), 
W Her (0.86), and T Her (0.15) which do not appear to have any type of Ca II IR triplet 
features. 

At this point, we call attention to three of the hydrogen Paschen lines wfiicfi lie close to 
eacfi Ca II IR triplet fine: 8502.4 A (4.4 A redward from the Ca II A8498.0 line), 8545.3 A 
(2.8 A redward of the A8542.1 line), and 8665.0 A (2.8 A from the A8662.1 line), the Pal3, 
Pal2, and PalO lines respectively. It has been shown by Gillct ct al. (1985) that even 
though the Pa-J line is in emission in the spectrum of the Mira variable S Car, the higher 
order Paschen fines near the Ca II IR triplet are neither seen in absorption nor emission, 
analogous to the Balmer fines. As such, it is unlikely that these higher order Paschen lines 
are affecting the Ca II lines in our spectra. 
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3.3. The TiO Bands and Ha Emission 

In addition to the Ca II IR triplet, we are interested in comparing the TiO 7 system 
molecular features to Ha emission as a function of phase. A qualitative comparison of Ha 
emission with the molecular features can be made for U Ori, R Leo, R CVn, R Ser, W Lyr, 
and X Cyg, stars with phase dependent spectra that also show the Ha emission line in at 
least one spectrum. Weak Ha emission is seen in U Ori at phase 0.29 on 15 Feb 1998. 
However on the previous pulsation cycle Ha emission is not obvious at either phases 0.30 or 
0.01, although there may some very weak emission. Stronger Ha emission is seen in R Leo 
(phase 0.34), R CVn (0.09, 0.17, and 0.39), R Ser (0.26), W Lyr (0.94), and x Cyg (0.04). 
As mentioned earlier in the paper. Ha is notorious for being observed as a weak emission 
feature when H/5, H7, and H5 are strong. Surprisingly though. Ha emission was not seen in 
various spectra where we would expect to find it: R Tau (phase 0.16 and 0.99), R Leo (0.16 
and 0.25), V CVn (0.12, 0.27, and 0.06), R CVn (0.20), V Boo (0.12), and T Her (0.15). 

Merrill (1940) and Joy (1926) report that the TiO bands are regularly stronger at 
minimum than at maximum light in Mira variables. To measure this trend, we checked 
the variability of the TiO feature at 7054 A (bandhead) with respect to a portion of the 
flux uncompromised by TiO, VO, and ZrO. We integrated the flux of each spectrum in 
the 6995-7045 A wavelength band (non-TiO) and the 7060-7110 A band (TiO), each 50 A 
wide. We then divided the integrated flux of the TiO band by the non-TiO band. By doing 
this, any scattered light that may be in the spectra are effectively canceled out. Since our 
data set contains warm oxygen- rich (earlier than M6), cool oxygen-rich (M6-M8), and one 
MS Mira sampled sporadically over various phases, we only include the cool oxygen-rich 
Miras in Figure 2, which graphs the above mentioned flux ratio as a function of phase. The 
M6-M8 stars are selected here in order to minimize the spread of effective temperatures at 
maximum light which will influence the strength of the TiO bands. Ideally, one would want 
many observations of each star over a single pulsation cycle. However, with the limited 
sample we have, we feel that we can get an approximate test of TiO variation with respect 
to phase. There appears to be no obvious trend in variations in the TiO 7054 A band flux 
with respect to the non-TiO band flux. Due to this observation, variations seen in Ha as a 
function of phase must primarily result from variations in the Ha emission itself and not to 
variations in overlying TiO absorption. As a result, using Ha flux variations as a proxy to 
variations in the intrinsic He flux is valid from this analysis. 
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4. Discussion 

We wish to test the idea that the apparent lack of He emission at 3970 A when the 
other Balmer hnes are strong emission features is anticorrelated with the strength of the 
Ca II absorption line at 8662 A. As reported in the Introduction, this anticorrelation results 
from He photons being scattered by the Ca II H line out to the Ca II hne at 8662 A, causing 
this Ca II absorption hne to be filled in with respect to the other two Ca II IR triplet lines. 

We use Ha as a proxy for the He line. Following the same analysis presented by 
CL97, we determined a relative line strength, F, for Ha emission and the Ca II IR triplet 
absorption lines. Two points were selected on either side of the emission or absorption 
feature. The wavelengths of these points were kept constant for all measurements. The 
observed profile was integrated across the wavelength window defined by these two points 
resulting in an integrated flux fi. A straight line connected between these two points 
represent a pseudocontinuum and the integrated flux, /c, calculated for it. Then 

I-, fe — fc 

where F will be negative for absorption lines and positive for emission lines. The 
measurements were done for the stars in our sample where their spectra included 
wavelengths below 6563 A and above 8662 A; a total of 27 spectra. Figure 3 shows the 
relative line strengths of Ha, and the Ca II A8498, Ca II A8542, and Ca II A8662 absorption 
lines as a function of phase. The uncertainty in the measurements is about ±0.007. The 
relative line strength of Ha clearly shows large scatter near visible maximum, and is zero 
within the uncertainty from phase 0.5 to 0.8. This is consistent with Balmer emission lines 
becoming prominent near maximum visible light. The Ca II A8498 line strengths are zero, 
within the uncertainty of the measurements; variation is not observed in this line. The 
Ca II A8542 line does show some scatter near phase 0.0, and is zero after phase 0.2. The 
Ca II A8662 line shows more scatter than either of the other two Ca II IR triplet lines, 
particularly near phase 0.0. 

From Figure 3, it is difficult to see any correlation between the relative line strength 
of the Ha and Ca II IR triplet lines. However, we can plot the relative line strengths of 
Ca II IR triplet lines versus the strength of the Ha line to look for correlations. Figure 4 
shows the relative line strengths of Ca II A8498, Ca II A8542, and Ca II A8662 versus the 
relative hne strength of Ha. The data of each plot is hnearly fit and the results of the fits 
are drawn in the plots. The slopes of the Ca II A8498 and A8542 relative hne strengths 
versus Ha are 0.08 and 0.05, respectively. The slope of the linear fit of the Ca II A8662 
relative line strength versus Ha relative line strength, on the other hand, is 0.32, which is 
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significantly different than the other two Ca II IR triplet hnear fits. Furthermore, the Ca II 
A8662 versus Ho; relative line strength slope is positive, so that as the strength of the Ha 
line increases, the strength of the Ca II A8662 line decreases (i.e. becomes more positive) 
and even goes into emission. This is the effect we expect if He photons are being scattered 
by the Ca II H line out to the Ca II line at 8662 A, causing the Ca II A8662 absorption line 
to be filled in with respect to the other two Ca II IR triplet lines. 



5. Conclusion 

The 6200 A to 9000 A spectra of Mira variables taken at different phases support a 
possible anticorrelation between Ha emission and Ca II A8662 absorption as first suggested 
by CL97. Assuming that the He line strength variations are in phase with Ha, then the 
apparent anticorrelation between the strength of the Ha emission line and the strength of 
the Ca II A8662 line suggests that a fluorescence is taking place in the Ca II A8662 line with 
He serving as the pump through the Ca II H line. This type of fluorescence is common in 
Mira type variables. The strong Fe I (42) lines at 4202 A and 4308 A seen in Miras are well 
known fluoresced features; in this case, the ultraviolet Mg II h & k lines serve as the pump 
via an Fe I (UV3) transition (e.g., Bidelman & Herbig 1958; Willson 1976; Luttermoser 
1996). 

The next phase of this research program is to systematically determine T^q and log g 
as a function of phase for the Miras in our sample. A detailed description of the LTE model 
synthetic spectra is given by Piontek & Luttermoser (1999). 
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Fig. 1. — Spectra of the 20 Mira variables. The name of the star, date, and phase, are given 
on each spectrum. The flux is normalized to one. Above each set of spectra are markers for 
the major spectral features, and terrestrial oxygen. Ha is weak, or not seen in most of the 
spectra. Markers enclose regions of individual spectra where no data was taken. Spikes due 
to cosmic ray hits have not been removed. 

Fig. 2. — The integrated flux ratio of the 7060-7110 A wavelength band (A) to that of the 
6995-7045 A band (B) plotted as a function of phase. Band A contains TiO opacity whereas 
band B is free of TiO, VO, and ZrO bandheads. This plot contains only data from the M6, 
M7, and M8 spectral-type Miras. No apparent trends are seen in the data, which indicates 
that variations in the Ha flux result primarily from intrinsic flux variations in Ha and not 
from varying overlying TiO absorption. 

Fig. 3. — Relative line strengths of Ha (flllcd square), Ca II lines at 8498 A (circle), 8542 A 
(triangle), and 8662 A (diamond), as a function of light- variation phase. Variation in the 
Ca II lines at 8498 A and 8542 A mimic each other, whereas the Ca II line at 8662 A does 
not follow the same trend. The uncertainty of the measurements is ±0.007. 

Fig. 4.— The relative hue strengths of Ca II A8498, Ca II A8542, and Ca II A8662 versus 
the relative line strength of Ha. The data of each plot is linearly flt and the results of the 
flts are drawn in the plots. Only the Ca II A8662 versus Ha plot shows a signiflcant slope 
which implies a correlation between the occurrence of the two features. 
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